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by 

Theodore L. Felsentreger 


SUMMARY 

A classification is presented of lunar satellite orbits perturbed by the earth 
and by the non- sphericity of the lunar gravitational field (the second zonal har- 
monic only). The orbits are arranged in three categories — those in which the 
argument of perilune circulates, those in which it "librates" about some odd 
multiple value otrr/2, and a class of "transition" orbits which belong to neither 
of the other two classes. 

The discussion is kept as general as possible by holding the semi-major axis 
arbitrary. In addition, two examples are presented to clarify the analysis. 
Graphs for these cases show the permissible ranges of the two parameters in- 
volved (integrals of the satellite motion) for the orbit classes. 



CLASSIFICATION OF LUNAR SATELLITE ORBITS 


INTRODUCTION 

Several studies have been made of the long-term effects on artificial satel- 
lite orbits disturbed by a third body (references 1, 2, 3 and 4). In addition, 

Lidov (4) has extended the analysis to include effects caused by the non- sphericity 
of the gravitational field of the central body (the second zonal harmonic only) , in 
the particular case when the inclination of the disturbing body's orbit plane to 
the equatorial plane of the central body is 90°. However, this simplification is 
not applicable when the central body is the moon, the angle in question being 
about 6°41*. 

The orbital motion of an artificial lunar satellite has already been studied 
by several people, among whom are Kozai (5) and Giacaglia, et al. (6). However, 
adequate analytical representations for the long-period and secular effects 
caused by the earth have proved particularly troublesome to obtain — for in- 
stance, these two papers present solutions involving elliptic integrals. Later 
solutions by Frost (7) and Fisher and Felsentreger (8) (which do not involve 
elliptic integrals) depend upon whether the motion of the argument of perilune 
is circulatory (secular) or libratory (periodic). It was therefore deemed ap- 
propriate to define the regions of circulation and libration. 

The disturbing forces which will be considered here, then, give rise to 
long-period and secular effects caused by the earth and secular perturbations 
due to the second zonal harmonic of the lunar gravitational field. 

Orbits will be divided essentially into two classes — those in which the 
argument of perilune either circulates or librates. The discussion will be limited 
to a presentation of the ranges of permissible values for two parameters (integral 
constants), which of course depend upon initial conditions. Values of the param- 
eters for orbits which fall into neither category (’’transition" orbits) lie on the 
boundaries of the ranges. 

A list of symbols appears in Appendix A. 
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EQUATIONS OF MOTION AND CONSTANTS OF THE MOTION 


The equations of motion for 77 = /l - e2 and the argument of perilune g 
are (6) 


= _Jj? — l n (1 - V 2 ) (1 - cos 2 i) sin 2g 
dt 4 q 


( 1 ) 


dg _ 
dt 


— — - 2 r ] 2 + 5(i7 2 -cos 2 i) sin 2 g 

2 q 77 


3 n . , c 2 • ^ 

— a. (1 - 5cos 1) . 

4^4 

77 


The two integral constants of the motion are readily found to be 


a = 7] 2 cos 2 i 


( 2 ) 


c - (1 - V 2 ) (1 -~sin 2 isin 2 g)-“A 


1 . 1 - 3 cos 2 i 


where 


A= 



In the case where the moon’s second zonal harmonic J 2 is neglected (i.e., 

A = 0) , the parameters become the a and /3 of Lorell ( 1 ), who has presented a 
classification of orbits disturbed only by a third body. However, the analysis 
is considerably more complicated when J 2 is not considered negligible. 

To begin the study, it will be assumed that for all orbits the argument of 
perilune must reach some odd multiple of n/ 2 . Therefore, the maximum c, a 
region compatible with elliptic orbits is that for which sin 2 g = 1 is allowable. 
If the argument of perilune is to circulate, then g must eventually become 0 . 
Hence, the region in the c , a plane describing "circulatory" orbits is that 
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for which both sin 2 g = 0 and sin 2 g = 1 are allowable. The non-intersection 
of the two regions gives "libratory" orbits. 

VALUES OF a AND c FOR WHICH SIN 2 g = 1 

For sin 2 g = 1, c can be rewritten as 

3ii 7 - (3 + 5a) r ) 5 + 5ari 3 - i At ) 2 + Aa 

c = : : 3 (3) 

2 r; 5 

which, for ~o = l (i. e ., e = 0), describes the line 


c 


_ I A (1 - 3a) . 


(4) 


For a particular value of rj , a achieves its greatest value when cos 2 i = 1. Re- 
placing Tj 2 in Equation (3) by a, one obtains 


c 



(5) 


which can be rewritten as 


c = - — A(l-3a) -- 
6 


0 1/2 -l) ^l + lA)a 2 +^l+lA^ 


|a 3/2 + i A a + i A a 1/2 + — A 
3 3 3 j 


, 3/2 


( 6 ) 


t 


The last term in Equation (6) is always £0, indicating that Equation (5) describes 
a curve in the c, a plane to the right of the line c = -(1/6) A(1 - 3 a), and which 
asymptotically approaches the axis a = 0. Equations (4) and (5) intersect at 
a = 1|C =(l/3)A. Therefore, the region to the right of the line of Equation (4) 
and bounded by Equation (5) and a = 0, gives permissible values of a and c . 


There remains the region to the left of the line to be explored. Equation 
(3) can be rewritten 


7 2 (2c + 3) T 7 3 - 3r) 5 + -j A 
5 (t ? 3 - 77 s ) + A 


( 7 ) 


3 


holding c constant and taking the derivative with respect to V yields 


2 77 1 1 577 10 -30 t 7 8 + 5(2c + 3) 77 s - 8A77 5 + 5A (c +^-) t} 3 +-^A 

da L 

d7 ?~ [5(77 3 - 77 s ) + A] 2 

It is desirable to describe the roots of 



f (77) = 1 577 10 


3O77 8 + 5(2c + 3 ) T7 6 - 8A77 5 + 5A 



+— A 2 


(9) 


in order to determine possible maximum values of a. 

Now, a = 0 at 77 = 0 . Since f( 0 )=(l/ 3 )A 2 > 0 , a increases from 0 until 77 
assumes the value of a root of f (77) . Suppose, for the moment, that Vi is a root 
of f (77) . If it can be shown that a ( 1 ) < a (tj 1 ) , then one can conclude that a 
attains a maximum value either at 77 x or at some other root of f (77) in 0<77<1. 

From Equation ( 9 ), then, 

- 15 77™ + 30 77® - 1577® + 8 A 77® - y A773 _ I A 2 

c = : 

5 hi (2 77 3 + A) 

Substitution of Equation ( 10 ) into Equation ( 7 ) yields (for 77 = 77^ 

77 2 [- 30 77™ + 30 77* + Art* + 5A77® + a 2 ] 
a(7 ?1 ) = Jl 1 1 1 , 

5 [5(77® - 7 ?®) + A] (2773 + A) 


( 10 ) 


( 11 ) 


and, for 77= 1, 


a(l) 


3077™+ 6O77 8 ( 9 -A) t }\ + 16 A 77® _ |- A (8 -A) 77^-^ A 2 
5 Atj 3 (2773 + A) 


( 12 ) 
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Then, 


a(l) - a(7 ?1 ) = 


0?! - l) 2 g(T7i) 

> 

15A7; 3 (2tj 3 + A) [5(7^ - 77 ®) + A] 


(13) 


where 


g(7 ?1 ) = 90(5 + A) 7?| 3 + 180(5 + A) 77 j 2 + I 8 OA 77* 1 - 180 (5 - A) 77 J 0 - 10 (45- 13A)t7® 
-A (250 + 3A)t7®-2A(290+3A)77j - A(290 + 49 A) 77 ® - 92A 2 77 * - lOOA^ 


-A 2 (50 + 3A)77 3 - 6 A 3 77 J - 4A 3 t] 1 - 2A 3 . 


(14) 


The Theorem of Vincent and Descartes’ Rule of Signs establishes that g (V !> 
has no root in 0 < 77 x < 1 (see Appendix B). Since g (0) = - 2A 3 < Oand g(l) = 
-15A(A 2 + 20A + 24) < 0, one may conclude that 


a (1) - 0(77^ < 0 


for 0 < 77 x < l; equality occurs only when r n l = 1. 

It must now be shown that a ( 77 ^ <1 for any root V 1 of f (v). From Equa- 
tion ( 11 ), 


a ( 7 7 1 ) 1 1 <" •> h(77 1 ) = -3077 j 2 + 3O77J 0 + 5O77® + A77 J - 50t7j 


+ 30 A 77 3 - 35 At 7 3 + A 2 7) 2 - 5 A 2 < 0. 


Again, it can be shown that h ( 77 ^ has no root in 0 1 Vi 1 1 (see Appendix B). 
Since h(0) = - 5 A 2 < 0, it has been established that 0 .( 77 ^ < 1 for any root o f f ( 77 ). 
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Finally, it must be demonstrated that, for any value of c less than 
-(1/6) A (l-3a), f (17) has at least one root. From Equation (9), f (0) = (1/3) A 2 > 0. 
Existence of a root between 0 and 1 is assured if f (1) < 0 - hence 


f(l) < 0 < 


;>c <A_( 4 -A) 
15(2 + A) 


However, when c = A( 4 -A )/15 (2 + A), 17 = 1 is a root of f (77), and the point a = 
(6+A)/5(2+A),c = A( 4 -A )/15 (2 +A) lies on the line c = -(l/6)A(l - 3a). Thus, 
c < A (4 - A)/l 5 (2 + A)=> f (77) has at least one root 77 x in 0 < r) l < 1 . in addi- 
tion, since a (1) < a (77 x ) for 77 x / 1, the point a (77 T ) , c lies above the line c = 
-( 1 / 6 ) A (1 - 3 a). 


Hence, the curve whose parametric equations are (for 0 < ??i £ 1) 

- 1577 j° + 3O778- 15776 + 8 A 775 _^At 73 _Ia 2 

c(7?i) = ^ 

5 r,l (2 + A) 

17 2 (- 3O77 j° + 3O77® + A77® + 5A77^ + A 2 ) 

°<Vi) = 

5 [ 5 (t 7 | - tjI) + A] ( 2 t 7 f + A) 


(15) 


represents the upper boundary to the permissible c, a region lying to the left 
of the line c= -(1/6) A (1 - 3 a). It intersects this line at the point c = 

A(4 - A)/15 (2 +A), a = (6 +A)/5(2 + A), and approaches the a = 0 axis asymptotically. 

VALUES OF a AND c FOR WHICH SIN 2 g = 0 


For sin 2 g = 0, Equation (2) becomes 


c 


- V 



+ — Aa 
2 


(16) 


which, for 77 = 1, again becomes the line c =-(l/6)A(l - 3a). in addition, for 
cos 2 i= 1, Equation (16) becomes 

Ia 

c = 1 - a +— , 
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which describes the same curve as for sin 2 g = 1. Therefore, the permissible 
region to the right of c = -( 1 / 6 ) A (1 - 3a) is the same as for sin 2 g = 1 . 

Now, Equation (16) can be rewritten as 


a = | V 2 + ( c - 1 ) 7 ? 3 + -jr aJ > 
whose derivative with respect to 17 (for constant c) is 

^ = |,[7V*5(c-1)V+}a] 

Thus, an analysis of the roots of 

F(7j) = 777 s + 5 (c - 1)77 3 + I A 


(17) 


(18) 


(19) 


is in order. 

Obviously, since F( 0 ) = (l/3)A>0, a increases from a value of 0 at '0 = 0 
until a root of F (77) is reached, assuring the existence of a relative maximum 
for a in 0 < 77 < 1 (providing F (77) has a root). Hence, assuming 77 1 to be a 
root of F (77), it again behooves one to show that a (1) < a 

When v = V 1 is a root of F (77), Equation (19) gives 


c 



( 20 ) 


Equation (17) then becomes (for 77 = 77 x ) 


a (-7 x ) 


V 1 (- 4vf + A) 


5A 


( 21 ) 
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for 17 = 1, 


-42 vf + 5(6 + A)vf - 2 A 

a (l) = : • (22) 

15 A V f 


Then, 


a(l) - a(7] 1 ) 


(Vi - I ) 2 GQh) 

15 A 77 j 


(23) 


where 


G(t 7 j) = 1 277 ® + 24t) J + 36i7j+48i7®+60t7^ + 3 (10 - A) 77 ^ - 6A77j -4At7j -2A. (24) 


G (tjj) has exactly one positive root. For A < 14, this root is in 0 < 77 x < 1; for 
A = 14, the root is 77 x =1; and for A > 14, G (77 x ) has no roots in 0 < Vi 1 1 
(see Appendix B). Therefore, a (1) - a (77^ < 0 for A > 14; for A < 14, designat- 
ing the root of G (77 x ) by 77 * , 


0.(1) - a(77 1 ) 


(V t - l ) 2 (^i " V\) G ‘(^i) 
15 A 77 ^ 


where G* (77^ > 0. Thus, a (1) - a (77^ < 0 when 77 j < 77* - equality occurs (for 
A < 14) only when 77 j = 77 * . 

The foregoing has also shown that the point c , a (77 t ) lies above the line 
c = -(1/6) A (l-3a) when 0 < V j < 1 (for A > 14) , and when 0 < 77 x < 77 * (for A < 14) . 
For the two cases respectively, the points 


c = - 




-4+A , 

a = and c = 

5A 


-7 17 * s + 5 t 7* 3 — jA 


577; 


*3 


t ?; 2 


(-477; 


*5 


a = • 


5A 


+ A) 


are on the line . 
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As before, now, it must be shown that a. (7^) < 1 for any root 77 x of F(77) . 
From Equation (21), 


a(?? 1 ) < 1 <s=> H (i7j) = -4 77J + A 77 2 - 5A 1 0. 


But H (T? 1 )has no root in 0 < Vi 1 1 (see Appendix B). Since H (0) = -5 A < 0, 
a(->7 1 ) < 1 for any root 77 x of F(77). 

Finally, in the case where A > 14, 


c < 



=0 F(l) = 5c + 2 + Aa<0. 


Also, F (0) = (1/3) A >0, so F (77) has at least one root in 0 < 77 <1. The existence 
of at least one root ofF (77) in the case where A < 14 will not be proved here — 
however, it can be proved after choosing a value of A < 14. 

To summarize, then, the curve whose parametric equations are 




a (Vl) = 


-lrj\ + 5?7^ A 
577? 

V] (- 4 i 7 ? + A ) 

5A 


(25) 


represents the upper boundary (for sin 2 g = 0) to the permissible c , a region 
lying to the left of the line c = -(1/6)A (1 -3a). When A > 14, the parameter V x 
may take on all values in 0 < 77 1 <1. However, when A < 14, 77 x is restricted to 
the range 0<77 1 <77*<1, where V 1 is the root of G (r^) defined by Equation (24). 

Figure 1 depicts the c , a regions defined in the previous two sections. 

k 
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EXAMPLES 


The two cases to be presented here correspond to values for the semi- 
major axis of the satellite’s orbit of 2 and 7.4822577 moon radii. Then, the 
values for A are 164.97081 and .22510948, respectively. For the first case, 
Equations (4), (5), (15), and (25) become 


c = - 27.495135 (1 - 3a) 


(4’) 


, , 54.990270 

c = 1 - a + 775 


( 5 ’) 


ft ' 


10 


-1 5T? j° + 30r)i - 1577® + 1319.7665 t)5 _1 099.8054 771 - 9071.7893 

5vf (2tj j + 164.97081) 


( 15 *) 




c (Vi) 


( - 30 77 + 30t)® + 164.970817)5 + 8 24.85405 77 * + 27215.368) 

5 [5(7)3 _ ^5) + 164.97081] ( 2rj\ + 164.97081) 

c , (0 < Vi— 1) 

— 7 77 5 +5?)j _ 54.990270 


5t) 


3 

1 


(25') 


Vi (-47)5 + 164.97081) 

C (7) ) = 

824.85405 (0 < Vl < 1 ) 

For A = .22510948, the root of G (t) 1 ) is 77 * = .25110445; therefore, the equations 
defining the boundaries of the circulation and libration regions are 


c = - .037518247 (1 - 3a) 


(4") 


c = 1 - a + 


.075036494 

n 3/2 


(5”) 


cOli) = 


- 15 77 j° + 30 77 ® - 15t)® + 1.800875847)5 _ 1.50072988 - .01689143 


5t)5 ( 2 77 j + .22510948) 


(15") 


a 0 ?i> 


7)2 (- 30 7 ) 1 ° + 307)8 + .22510948 v\ + 1. 12554740 77 * + .05067428) 
5 [5(7)3 _ 7)5) + .22510948] (2t)3 + .22510948) (Q < < 1} 


C (.V X ) = 


- 77)5 + 5 7)5 + .075036494 
57)5 


(25") 


a = 


7)2 (- 4 T)® + .22510948) 
1.1255474 


(0 < 77 , < .25110445) 
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Tables 1, 2, and 3 (Appendix C) list values of c and a for points on the 
curves defined by the preceding equations. Figures 2, 3, 4, and 5 (Appendix D) 
show the circulation and libration regions for the two cases. It is readily ap- 
parent that the libration region for a =2 1^ is considerably smaller than that 
for the case a = 7.4822577 R c - this is a consequence of the fact that the part 
of the disturbing function dependent upon J 2 (which is purely secular) is more 
significant for smaller values of a . 


CONCLUSIONS 

Given a set of initial conditions for an artificial lunar satellite, it can be 
determined from a graph similar to Figure 1 whether the orbit falls into the 
circulation case or libration case. If circulation, the method outlined by 
Frost (7) can be used to compute the long-period and secular effects caused by 
the earth and the secular effects due to the moon’s second zonal harmonic. On 
the other hand, if the orbit is libratory, the method presented by Fisher and 
Felsentreger (8) may be employed. 
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APPENDIX A 


SYMBOLS 

a= semi-major axis of satellite’s orbit 
e= eccentricity of satellite's orbit 

i = inclination of satellite's orbital plane to moon's equatorial plane 
g= argument of perilune of satellite's orbit 
n= mean motion of satellite 
t= time 
r}= V 1 - e 2 

R c = mean radius of moon - 1738 km. 

a c = semi-major axis of moon's orbit - 221.17376 R^. 

J 2 = second zonal harmonic coefficient of moon's gravitational potential 
— 2.41 x 10 4 

1 t mass of moon ^ t A100 
q = 1 + — l.OIZo 

mass of earth 
a i = (V2) (a/a c ) 3 

a 2= K J 2 / a2 

A=-^- 

a, . 

l/q 

a, c = constants of the motion 
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APPENDIX B 


1. Proof that g('n 1 ) has no root in 0 < 17 < 1 
From Equation (14), 

g (77 1 ) = 90 (5 + A) 7i\ 3 + 180(5 + A) 77 j 2 + 180 A 77} 1 - 180(5 -A) 77 - 10(45- 13A) 77 * 

- A(250 + 3 A)? 7 j -2A(290 + 3A) 77 7 -A(290 + 49A) 77 ® -92A 2 t ? 3 -100A 2 77 * 

- A 2 (50 + 3A) 77 3 - 6A 3 t ] 2 - 4A 3 t) x - 2A 3 , 


which, by Descartes’ Rule of Signs, has either one positive root, or possibly 
three (in the case 45/13 < A < 5) . Replacing v 1 by l/l + X (X > 0) transforms the 
equation into one having exactly the same number of positive roots as those of 
g (t?j) in 0 < < 1 . Thus, solving g ( 17 x ) = 0 for 0 < 77 x < 1 is equivalent to 

solving the following equation forX >0 (the coefficients of X 13 , X 12 f . . . X s are 
all obviously negative, so there is no need to write them down explicitly): 


0 = a 13 X 13 + a 12 X 12 + a n X 11 + a 10 X 10 + a g X 9 + a 8 X 8 + a 7 X 7 + a 6 X 6 + a 5 X s 

-10(45 + 1997A+ 3136A 2 +602A 3 )X 4 -3(900 + 7850A + 7 1 39A 2 + 934A 3 )X 3 
-5(1080 + 31 58A+ 191 5A 2 +177A 3 )X 2 - 10(360 + 516A + 253A 2 + 17A 3 )X 
- 15A(24 + 20A + A 2 ), 

where a x 3 , a 1 2 » • • *> a s < 0. This equation obviously can have no positive 
solutions - hence, g (Vi) has no roots in 0 < 77 x <1, 

2. Proof that h ( 77 x ) has no root in 0 < V\ 1 1 


h( 7 7 1 )= - 30 77 ^+ 3077^° + 5 O 77 8 + A 77 J - 507? 8 + 30A77j - 35krj 3 + A 2 77 2 - 5 a 2 . 
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Since h (0) = -5 A 2 < 0 and h(l) = -4 A(l+A)< 0, it remains to be shown that 
h(?] t ) has no roots in 0 < t) x < 1. Equivalently, it must be shown that the trans- 
formed equation, after replacing 77 x by l/l +X (X>0), has no positive roots. The 
transformed equation is 


0 = -5 A 2 X 12 - 60A 2 X 11 - 329A 2 X 10 - 5A(7 + 218A) X 9 - 45A (7 + 54A) X 8 

-30 A(41 + 128 A) X 7 - 10(5 + 273A + 441A 2 ) X 6 - (300 + 3779A + 3708A 2 ) X s 
-5(140 +671A + 453A 2 ) X 4 -20(40 + 94A + 49A 2 ) X 3 -5(84 + 124A+ 57A 2 ) X 2 
-10(4 + 10A+ 5A 2 ) X - 4A (1 + A), 


which has no non-negative solutions. Consequently ,h (77 x ) has no roots in 
0 < 77 x <1. 

3. Roots of G (17 in 0 < -rq < 1 


G(t 7 x ) = 1 2 t 7 8 + 2477 7 + 36 r,\ + 48t7 3 + 6O77 4 + 3(10 - A)t 7 3 - 6A77 2 - 4 At 7 x - 2A. 


Finding the roots of G (77 x ) in 0 < 77 x < 1 is equivalent to finding the positive 
roots of 

b g X 8 + b 7 X 7 + b 6 X 6 + 5(1-47A) X s + 35(6 -llA)X 4 + 6(98 -67 A) X 3 
+ 20(42 - 1 3A) X 2 + 5(126 - 19A) X + 15(14 - A), 


where b g , b ?> b 6 <0 and 77 x = 1/1 + X. Descartes' Rule of Signs establishes 
that the polynomial in X has exactly one positive root for A < 14, and no positive 
roots when A > 14 . Consequently, G (77^ has exactly one root in 0 < Vi 1 1 for 
A 1 14 , and no roots in 0 < 77 x < 1 when A > 14 . Since G(l) = 15 (14 - A), the root 
is 1 when A = 14 , and is in 0 < 77 x < 1 when A < 14 . 
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4. Proof that H(t? 1 ) has no root in 0 < \ < 1 


H-iV-i) — 4t)j + A17 2 - 5 A. 


H(l) = - 5A < 0, so 1 is not a root. 

Setting Vi = 1/1 +X , the polynomial set equal to zero is 


0 = -5 AX 7 - 3 5 AX 6 - 104 AX 5 -170AX 4 - 165AX 3 - 95AX 2 - 30AX-4(1 +A), 


which obviously has no positive roots. Therefore, H (t^ ) has no roots in 

0 < T7j < 1. 
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Tables 
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Table 1 


1 

3 A 

c, a values for c = 1 - a + 


a 

c (A = 164.97081) 

c (A = .22510948) 

. i 

54.990270 

.07503649 

.9409 

60.310987 

.14131618 

.8836 

66.323153 

.20674184 

.8281 

73.144811 

.27147455 

.7744 

80.918967 

.33570943 

.7225 

89.819971 

.39968440 

.6724 

100.061843 

.46369149 

.6241 

111.909208 

.52809180 

.5776 

125.619822 

.59333530 

.5329 

141.824080 

.65998744 

.4900 

160.831490 

.72876529 

.4489 

183.386990 

.80058712 

.4096 

210.361630 

.87664150 

.3721 

242.896060 

.95848490 

.3364 

282.503110 

1.04818164 

.3025 

331.217530 

1.14850822 

.2704 

391.818790 

1.26325735 

.2401 

468.169490 

1.39769967 

.2116 

565.741440 

1.55930176 

.1849 


1.75887216 

.1600 


2.01244520 

.1369 


2.34448300 

.1156 


2.79353121 

.0961 


3.42266381 

.0784 


4.33980761 
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Table 2 


— 1 5t?J° + 30^7® - 1 5t7| + 8At^ _ yA7]3 - -A 2 

c, a values for c = — - ~ 

5t?3 (2^3 + A) 



n — 

Tj\ (-307,*° + 3 O 17 ® + 

Kt}\ + 5At,3 + A 2 ) 



(X — * 

5 [5(7,3 _ fjS) + A ] (27,3 + A ) 



c(A= 164.97081) a(A = 164-9708 1) 

c(A= .22510948) 

a(A= .22510948) 

1 

- 10.60285 

.20479126 

.02545983 

.55953289 

.95 

- 12.58625 

.18367060 

- .00179364 

.45293847 

.90 

- 14.99181 

.16402945 

- .05465415 

.36267636 

.85 

- 17.95302 

.14574676 

- .12887024 

.28707338 

.80 

- 21.65669 

.12872644 

- .22039516 

.22453776 

.75 

- 26.36940 

.11289356 

- .32639709 

.17250927 

.70 

- 32.48019 

.09819074 

- .44028987 

.13267987 

.65 

- 40.57145 

.08457501 

- .56180739 

.10053275 

.60 

- 51.54089 

.07201504 

- .68716695 

.07578677 

.55 

- 66.82003 

.06048873 

- .81440876 

.05716197 

.50 

- 88.78450 

.04998109 

- .94307939 

.04342257 

.45 

- 121.56804 

.04048246 

-1.07559004 

.03340831 

.40 

- 172.78868 

t 

.03198710 

-1.21995835 

.02596810 

.35 

- 257.51846 

.02449195 

-1.39568007 

.02022638 

.30 

- 408.39139 

.01799570 

-1.64784077 

.01545777 

.25 

- 704.97551 

.01249808 

-2.08708494 

.01126196 

.20 

-1375.89282 

.00799932 

-3.02826351 

.00753283 
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Table 3 




-7v[ 

+ Srjf - 


v\ (~4v\ + A) 



c, a value 

s for c - 



a - 





5't]3 


5A 


ZLi 

c (A =164.97081) a (A =164.97081) 



c (A = .22510948) 

a (A =.22510948) 

1 

- 11.39805 

.19515066 

.25110445 

- .03612413 

.012386983 

95 

- 13.09109 

.17711353 

.240 


- .16623741 

.011357005 

90 

- 15.22049 

.15968058 

.230 


- .30750281 

.010458998 

85 

- 17.91999 

.14294541 

.220 


- .47716070 

.009591355 

80 

- 21.37657 

.12698302 

.210 


- .68222360 

.008755992 

75 

- 25.85696 

.11185269 

.200 


- .93191234 

.007954511 

70 

- 31.75030 

.09760064 

.190 


-1.23851181 

.007188233 

65 

- 39.63903 

.08426227 

.180 


-1.61862795 

.006458243 

60 

- 50.42092 

.07186425 

.175 


-1.84307067 

.006107137 

55 

- 65.52751 

.06042617 

.170 


-2.09506996 

.005765417 

50 

- 87.33443 

.04996211 

.165 


-2.37891665 

.005433167 

45 

- 119.97545 

.04048188 

.160 


-2.69973129 

.005110460 

40 

- 171.06859 

.03199205 

.155 


-3.06365715 

.004797361 

35 

- 255.68588 

.02449688 

.150 


-3.47810698 

.004493928 

30 

- 406.46133 

.01799894 

.145 


-3.95208005 

.004200211 

25 

- 702.96296 

.01249970 

.140 


-4.49657219 

.003916254 

20 

-1373.81275 

.00799994 

.137 


-4.86261985 

.003750580 




.133 


-5.40369117 

.003535184 




.130 


-5.85447414 

.003377770 




.128 


-6.17897523 

.003274799 




.127 


-6.34899307 

.003223906 




.125 


-6.70561192 

.003123305 




.123 


-7.08586163 

.003024286 
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APPENDIX D 
Graphs 
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-Libration and circulation regions for A = 164.97081 





-1 0.60285 ,.20479 1 26) 



Fig. 3-Libration region for A = 164.97081 





ig. 4-Libration and circulation regions for A = .22510948. 





Fig. 5-Expanded view of curve N # from Figure 4 (A = .22510948). 




